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thin-layer  dosimeters  prepared  either  as  self-supporting  films  of  powder  in 
a  polyimide  matrix,  or  on  substrates  of  LiF  single  crystals  or  borosilicate 
glass.  It  is  demonstrated  that  all  relevant  opcical  and  thermal  properties 
of  the  dosimeters  can  be  determined  by  these  methods  and  that,  based  on  this 
knowledge,  the  expected  thermoluminescence  response  of  a  given  configuration 
can  be  simulated  as  a  function  of  a  number  of  experimental  parameters. 

For  the  investigation  of  the  dependence  of  the  thermoluminescence 
emission  yield  for  a  given  dose  on  the  heating  rate  and  for  experimental 
tests  of  the  developed  theory  of  heating  continuous  dosimeter  configurations 
with  a  beam  of  uniform  intensity  profile,  a  new  thermoluminescence  laser 
heating  facility  was  designed  and  constructed.  It  features  an  advanced 
closed-loop  laser  control  electroncis,  a  40  Watt  laser,  and  a  uniform  round 
beam  of  3  mm  diameter  or  a  square  one  (selectable)  of  3  x  3  mm2  area.  This 
apparatus  is  described  in  Appendix  C. 
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The  purpose  of  this  research  effort  is  the  investigation  of  the  basic 
characteristics  of  two  differnt  types  of  laser-stimulated  luminescence 
phenomena  in  view  of  potential  applications  in  the  dosimetry  of  ionizing 
radiation  and  of  radiation  imaging. 

The  first  of  these,  laser-heated  thermoluminescence,  was  the  subject  of 
the  research  described  in  the  present  Annual  Report.  The  second 
phenomenon  will  be  the  main  subject  of  our  studies  in  the  next  contract 
period. 

The  heat  transfer  mechnanism  from  a  CO2  laser  beam,  having  a  Gaussian 
intensity  profile,  to  a  semi-infinite  thin  thermoluminescence  phosphor 
layer  and  the  resulting  luminescence  response  has  been  studied  in  detail. 
Experimental  and  computational  methods  were  developed  for  the 
chacterisation  of  a  given  thermoluminescence  dosimeter  configuration.  All 
its  relevant  optical  and  thermal  properties  can  be  determined  by  these 
methods.  Based  on  this  knowledge,  the  expected  thermoluminescence 
response  can  be  computer-simulated  as  a  function  of  a  number  of 
experimental  parameters.  This  in  turn  provides  a  tool  to  optimize  the 
dosimeter  design  for  applications  in  dosimetry. 


I.  Introduction 

In  recent  years  the  feasibility  of  heating  thermoluminescent 
materials  with  infra-red  laser  beams  has  been  demonstrated  [l-H]. 
This  new  thermal  stimulation  method  has  received  particular  attention 
in  thermoluminescence  dosimetry  (TLD)  of  ionizing  radiation  because  it 
holds  promise  as  a  solution  to  a  number  of  problems  associated  with  the 
measurement  of  small  doses  of  nonpenetrating  radiation  such  as  low 
energy  beta  rays  as  well  as  knock-on  protons  produced  by  fast  neutrons 
in  hydrogenous  radiators  [ 5— 7 ) -  Even  the  imaging  of  spatial  dose 
distributions  of  Y-rays  and  x-rays  has  been  contemplated  [ 1 , 2 ]  and 
experimentally  demonstrated  [8]. 

The  successful  development  of  practical  laser-heated  TLD  readers 
and  thermoluminescence  imaging  devices,  requires  detailed  knowledge  of 
the  heat  transfer  mechanism  from  the  laser  beam  to  the 
thermoluminescent  sample  and  the  resulting  spatio-temporal  temperature 
distribution.  Then  one  may  compute  and  predict  the  thermoluminescence 
response  (glow  curve)  of  the  sample  dosemeter  for  a  given  laser  heating 
beam  [9].  The  effects  of  such  Interrelated  design  parameters  as  laser 
power  density,  time  of  laser  exposure,  laser  beam  diameter  and 
intensity  profile,  sample  configuration  and  composition,  and  the 
expected  incandescent  background  emission  can  be  computer-simulated  and 
the  resulting  glow  curve  compared  with  experimental  results. 

The  list  of  experimental  parameters  is  narrowed  down  to  the  rew 
that  presently  appear  important  for  practical  applications.  For 
example,  the  laser  power  of  cost-effective  modern  RF  or  DC  discharge 
excited  C02  waveguide  lasers  or  RF  excited  unstable  resonator  lasers  is 
limited  to  less  than  100  Watts.  Laser  beam  profiles  are  normally 
Gaussian  intensity  distributions  or  the  characteristic  annular  cross 
section  ("halo")  that  is  possible  with  an  unstable  resonator  cavity. 
Uniform  beam  profiles  are  highly  desirable,  but  all  known  methods  to 
produce  such  beams  appear  to  be  impractical  in  TLD  applications  because 


of  the  associated  power  loss  or  the  required  custom  fabrication  of 
expensive  special  optics  [ 1 0 ] - 

Suitable  dosemeter  configurations  are  thin  layers  of  the  TLD 
phosphor,  either  in  the  form  of  a  self-supporting  film  composed  of  a 
high  temperature  polymer  mixture  with  phosphor  particles,  or  applied  to 
a  thin  substrate.  Both  may  take  the  form  of  discrete  circular  spots 
whose  diameter  is  smaller  than  the  full-width  at  1 /e-power  of  the  laser 
beam  or  may  be  prepared  as  continuous  layers  whose  area  is  much  larger 
than  the  laser  beam  diameter.  However,  even  the  widely  used  square  TLD 
chips  (e.g.  the  configurations  supplied  by  Harshaw/Filtrol  [  1 1  ] )  or 
round  pellets  [ 1 2 ]  can  be  successfully  heated  with  a 
COj-laser  beam  [ 1 3 ] • 

Small  discrete  dosemeters  can  be  heated  quite  uniformly  with 
Gaussian  beams  provided  the  laser  photons  are  not  too  strongly  absorbed 
(uniform  heating  in  the  direction  of  laser  propagation).  Strongly 
absorbing  dosemeter  layers  (surface  heating)  must  be  thin  enough  that 
the  thermal  response  time  [  1 4 ]  of  the  sample  is  much  shorter  than  the 
total  heating  time  required  to  release  all  trapped  carriers.  A 
necessary  condition  for  a  spatially  uniform  temperature  rise  in  both 

cases  is,  of  course,  that  the  diameter  of  the  discrete  3pot  is  much 

smaller  than  the  diameter  of  the  laser  beam.  Obviously,  a  considerable 
fraction  of  the  laser  beam  energy  is  wasted  in  this  mode  of  operation 
which,  for  this  reason,  will  probably  be  considered  only  in  special 
situations . 

Continuous  dosemeter  layers  are  of  practical  interest  because  of 
the  ease  with  which  they  can  be  fabricated  and  because  of  their 

potential  usefulness  in  a  number  of  the  new  applications  mentioned 

above.  A  localized  Gaussian  beam  can  be  employed  to  rapidly  heat  a 
small  3 pot  on  the  TLD  layer,  yielding  a  characteristic  glow  curve  that 
is  the  result  of  the  non-uniform  spatio-temporal  temperature 
distribution  T(x,y,z,t)  produced  in  this  way.  This  type  of 
thermoluminescence  response  is  completely  different  from  the 
conventional  glow  curves  obtained  by  uniform  contact  heating  (an 
example  is  shown  in  Fig.  9  below).  It  is  nevertheless  very  useful  in 
dosimetric  and  imaging  applications  of  thermoluminescence  [ 1 5 ] . 
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In  this  paper  we  present  the  theory  of  the  thermoluminescence 
response  curves  generated  with  Gaussian  laser  beam  profiles  in  a  number 
of  dosemeter  conf igurations  that  are  of  practical  importance  in 
thermoluminescence  dosimetry.  The  theoretical  approach  is  verified  by 
comparison  with  experimental  results.  In  addition,  experimental 
methods  are  presented  for  the  determination  of  all  relevant  thermal 
and  optical  properties  of  the  dosemeter  materials.  Thus,  the 
thermoluminescence  response  of  a  given  type  of  dosemeter  is  completely 
characterized  and  the  effects  of  such  design  parameters  as  layer 
thickness,  laser  power  and  beam  size,  pre-annealing  temperature  and 
duration,  substrate  material  and  its  thickness  can  all  be  assessed  by 
computer  simulation.  The  special  case  of  a  semi-infinite  LiF  slab  has 
been  reported  previously  [9]. 

II.  THERMOLUMINESCENCE  KINETICS  FOR  NONUNIFORM  SPATIO-TEMPORAL 

TEMPERATURE  DISTRIBUTIONS 

1  .  GENERAL  REMARKS 

Calculations  of  thermoluminescence  glow  curves  on  the  basis  of 
simple  electron  kinetic  models  and  for  spatially  uniform  heating  pose 
no  principle  problems.  While  only  special  cases  yield  analytical 
solutions  for  the  thermoluminescence  intensity  as  a  function  of  time  or 
temperature,  numerical  solutions  can  readily  be  generated  for  the 
coupled  nonlinear  rate  equations  which  describe  a  given  trap  model 
[16,17].  The  extreme  heating  rates  possible  with  lasers  have  been 
shown  not  to  invalidate  the  electron-statistical  foundation 
(Shockley-Read  statistics  [l8,l9])  of  these  phenomenological 
theories  [4]. 

Usually  a  number  of  different  physically  plausible 
electron-kinetic  models  will  give  satisfactory  fits  of  experimental 


glow  curves.  This  general  lack  of  uniqueness  renders  trap-level 
spectroscopy  by  monitoring  thermally  simulated  electron-kinetic 
relaxation  phenomena  a  nontrivial  task  [16,17].  Spatially  nonuniform 
heating  significantly  Increases  the  required  computational  effort  and 
further  complicates  this  situation.  For  this  reason,  laser  heating 
techniques  are  explored  herein  for  the  sole  purpose  of  assessing  their 
utility  in  the  dosimetry  of  ionizing  radiation.  No  attempt  is  made  to 
verify  a  given  simple  trap  model  or  to  extract  new  methods  which  might 
facilitate  trap-level  spectroscopy.  Thus,  all  that  is  required  for  our 
purposes  is  one  reasonably  plausible  model  description  for  the 
thermoluminescence  kinetics  of  a  suitable  phosphor  heated  with  a 
laser . 

The  calculations  presented  here  are  based  on  simple 
Randall-Wilkins  first-order  kinetics  [ 1 6  ] .  A3  the  example  for  a 
typical  phosphor,  we  choose  the  most  widely  used  dosimetry  material 
LiF:Mg,Tl  (TLD-100)  [ 1 2 ] .  Experimental  fits  of  its  prominent  glow 
peaks  (costumarily  labelled  peaks  2  through  5)  have  been  reported  by  a 
number  of  authors,  most  recently  by  McKeever  [ 20 ]  and  by  Vana  and 
Ritzinger  [21 ].  The  fact  that  the  numerical  fitting  parameters, 
notably  the  trap  depths  and  the  frequency  factors,  differ  significantly 
in  these  papers  only  underlines  the  problems  associated  with  curve 
fitting  as  a  method  to  obtain  basic  physical  information  from 
thermoluminescence  experiments  alone.  We  have  arbitrarily  selected 
McKeever 's  data  for  our  investigations. 

2.  COMPUTATIONAL  APPROACH 

Having  selected  a  group  of  glow  peaks  in  TLD-100  and  a  suitable, 
albeit,  not  unique,  fit  for  modelling  these  peaks  under  various 
experimental  conditions,  the  main  task  is  the  calculation  of  the  time 
evolution  of  the  temperature  ■  ("heating  program" )  for  each  volume 
element  dxdydz  at  location  (x,y,z)  of  the  sample. 

A  rich  literature  exists  in  the  field  of  laser  heating  of 
materials.  For  example,  laser  annealing  of  semiconductors  has  become 
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an  Important  manufacturing  tool  in  the  electronics  industry  [22-27]. 
Laser  absorption  calorimetry  of  very  transparent  optical  materials  also 
requires  the  knowledge  of  the  time-dependent  temperature  distributions 
produced  in  the  sample  [28].  Important  insight  in  the  laser  damage 
properties  of  high-power  metallic  and  dielectric  mirrors  and  other 
optical  components  has  been  gained  by  studying  the  theory  of  laser 
heating  [29, 30].  Comprehensive  treatments  of  heat  conduction  in  solids 
and  associated  boundary  value  problems  are  the  texts  by  Carslaw  and 
Jaeger  [31]  and  Ozisik  [32].  However,  despite  an  extensive 
bibliography  on  laser  heating,  the  problem  of  heating 
thermoluminescence  dosemeters  with  a  Gaussian  beam  turned  out  to  be  a 
unique  new  case.  It  is  perhaps  most  closely  related  to  the  one  solved 
by  Bernal  [28]  who  considered  Gaussian  beam  heating  of  extremely  weakly 
absorbing  cylinders  and  semi-infinite  slabs.  The  laser  beam  is,  in 
first  approximation,  unattenuated  when  traversing  the  entire  thickness 
of  the  slab.  Unfortunately,  this  special  case  is  of  little  interest 
here.  Weakly  absorbing  thermoluminescence  dosemeter  layers  are  not 
suited  for  practical  applications  because  of  inefficient  use  of  costly 
beam  energy.  Strong  to  moderate  beam  attenuation  in  the  active  layer 
is  characteristic  for  laser  heating  in  thermoluminescence  dosimetry. 
However,  this  problem  can  be  solved  analytically  by  the  Green's 
function  method  under  the  assumption  that  the  layer  is  semi-infinite  in 
radial  direction  (local  heating  of  a  small  3pot  on  a  large  area)  and 
that  no  significant  heat  loss  occurs  due  to  convection  and  radiation 
[9].  Convection  can  easily  be  eliminated  by  evacuation  of  the  sample 
chamber  and  radiative  losses  are  indeed  small  considering  the  fact  that 
the  maximum  temperatures  of  interest  do  not  exceed  700K  and  the  heating 
times  are  limited  to  a  few  hundred  milliseconds.  For  these  assumptions 
to  be  valid,  the  rate  of  energy  absorption  from  the  laser  beam  must  be 
significantly  greater  than  that  of  radiative  energy  loss. 

The  resulting  spatio-temporal  temperature  distribution  is 
cyl indr ically  symmetric  around  the  beam  axis.  To  calculate  the  total 
thermoluminescence  emission  as  a  function  of  time  after  onset  of  the 
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laser  exposure,  the  time-dependent  contributions  from  volume  elements 
2irrArAz  are  summed  by  numerical  integration  over  the  radial  coordinate 
r  and  the  direction  of  laser  beam  propagation  z.  Because  these  heating 
rates  are  nonlinear  in  time,  analytic  solutions  of  the  electron  kinetic 
rate  equations  for  each  of  these  inf initesmial  volume  elements  and 
their  individual  unique  heating  rate  are  unavailable  even  for  simple 
first-order  kinetics.  Therefore,  the  rate  equations  must  be  solved 
numerically  on  the  basis  of  the  individual  temperature  evolution  of 
each  of  these  volume  elements  and  their  contribution  to  the  emission 
intensity  I(t)  must  be  summed  up. 

The  typical  situation  encountered  in  heating  a  small  area  of  a 
continuous  serai-infinite  layer  of  thickness  L  by  a  Gaussian  laser  beam 
of  intensity 


I(r,z)  -  IQ  exp(-r i/wi-yz)  (1) 

is  schematically  depicted  in  Fig.  1.  The  symbols  chosen  are  che  full 
width,  2w ,  at  l/e-peak  power  of  the  beam,  the  absorption  coefficient, 
u,  the  radial  distance,  r,  from  the  axis,  and  the  distance,  z,  from  the 
front  surface  upon  which  the  beam  is  incident.  If  the  beam  is  turned 
on  at  t-0,  the  temperature  distribution  T(r,z,t)  evolves  as  a 
consequence  of  energy  absorption  and  thermal  diffusion,  resulting  in  a 
spatially  resolved  time-dependent  thermoluminescence  emission  pattern 
l„.  (r,L,t)  similar  to  the  special  case  shown  in  Fig.  2.  These 

i  L * 

photographs  were  obtained  by  heating  a  highly  absorbing  1mm  thick  glass 
slab  whose  back  side  at  z«L  was  coated  with  a  50um  thick  layer  of 
ZnS:Cu  powder.  Here  the  thermoluminescence  emission,  1^  (r,z,t),  is 
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not  due  to  direct  laser  heating  of  the  phosphor  itself.  Instead,  the 
thermoluminescent  layer  is  heated  by  temperature  diffusion  through  the 
glass.  Similar  experiments  were  described  previously  [2].  After 

appropriate  calibration  such  a  thermoluminescence  layer  can  be  used  to 
measure  the  time  evolution  of  Isotherms  during  laser  heating  of 
materials.  For  example,  the  diameter  of  the  dark  center  spot  in  Fig. 
2b  -  2c  corresponds  to  the  438K  isotherm  because  the  thermoluminescence 
emission  of  ZnS:Cu,  measured  by  conventional  uniform  slow  heating  at  a 
constant  heating  rate  of  15-3K/S,  ceases  at  this  temperature.  Note 

that  the  average  heating  rate  in  the  laser  experiment  (Fig.  2)  is  about 
the  same  as  can  be  discerned  from  the  appearance  of  the  dark  center 
spot  at  about  10s.  Other  isotherms  may  be  selected  a3  well.  For 

example  the  peak  emission  (region  of  maximum  brightness  in  Fig.  2) 
corresponds  to  approximately  38IK.  Still  other  emission  intensities, 
measured  relative  to  the  peak  emission,  may  be  chosen.  Actually,  the 
entire  temperature  distribution  present  at  a  given  time  after  onset  of 
heating  can  be  determined  by  comparison  of  the  normalized  emission 
intensities  in  Fig.  2  with  those  of  the  conventional  glow  curve 
measured  at  approximately  the  same  heating  rate. 

In  dosimetry  applications  of  laser  heating  the  emission, 

integrated  over  the  total  sample  volume,  is  measured  with  a 
photomultiplier  tube  without  spatial  resolution: 

ITL(t)  -  //drdz  Iy(r,z,t)  (2) 

Here  I  is  the  thermoluminescence  emission  intensity  per  unit  volume. 

As  a  first  step  in  the  theory  of  this  TL  emission  we  calculate 
T'r,z,t)  by  solving  the  thermal  diffusion  equation  [9]. 

72T(r,z,t)  ♦  k_1g(r,z,t)  -  a-1  dT/dt  (3) 

for  the  following  boundary  conditions: 

a)  The  slab  is  infinite  in  radial  direction  and  has  a 
thickness  L. 

b)  No  heat  loss  occurs  at  the  boundaries  z-0  and  z-L, 

i.e.  3T/3z-0. 
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This  assumption  is  valid  in  vacuum  in  the  absence  of  radiative  heat 
loss . 

The  Initial  condition  is  T(r,z,t)-To,  i.e.  room  temperature  when 
the  laser  beam  is  turned  on  at  t-0.  Then  the  temperature  obtained  from 
Eq.  (3)  is  the  actual  temperature  increase  above  T  .  The  notation  used 
in  Eq.  (3)  is  as  follows: 

k  «  thermal  conductivity, 
a  -  k/pc  -  thermal  diffusivity, 
p  »  mass  density, 
c  =  specific  heat, 

g(r,z,t)  -  laser  power  per  unit  volume  present  at  position  (r,z) 
during  laser  exposure  (source  function). 

The  case  of  temperature  dependent  k  and  ct  is  discussed  in 
Appendix  A. 

Correcting  for  reflective  loss  on  the  front  surface  at  z-0,  the 
source  function  g  becomes 

g(r ,z,t )  -  p  (1-R)  IQexp(-r2/w2  -  yz)  for  t  S  0  ^ 

-  0  for  t  <  0, 


where  3  is  the  reflectivity. 

We  solve  Eq.  (3)  with  the  initial  and  boundary  conditions  stated 
above  by  the  Green's  function  technique  [32;.  The  general  solution  is 

t  L 

AT  »  T(r,z,t)  -  T  -  —  /dT  /r'dr*  /  dz'  x 

°  t-0  r'-O  z'-O 

[ G ( r , z , t ; r ’ ,z' g ( r ' , z ' , t ) ]  (5) 


where  G  is  the  following  Green's  function: 
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G(ptztt;p* ,2' , t)  -  /d8e  °8  (t~T^8J0(Sr)J  (8r') 

8*o 


[—■*■  —  I  e  ar|m^t  T^cos(nniz)  cosCr^z')] 
L  L  m»  1 


(6) 


Here  J  is  the  zero-order  Bessel  function  and  n  -mirL- 1 .  After 
o  m 

integration  over  r’.z',  and  8,  one  obtains  with 

e» 

P  -  /I (r ,0)2irrdr-ir  w2Iq  (total  power  incident  on  the  front  surface  of 
o 

the  thermoluminescent  layer  at  z-0)  the  temperature  distribution  [9]: 


AT 


(1-R)aP 

kLir 


t 


■  /  dr 

T-O 


e~r*/[4a(t-T)  +  w2] 
^a(t-x)  +  w2 


x 


m 

2  l 

m»1 


-an2(t-x) 
e  m 

1+(n  /u)1 
m 


(1-e  wLcos 


mu) 


003 ( nz ) } 
m 


(7a) 


Note  again  that  Eq.  (7a)  is  the  temperature  in  excess  of  the  initial 

temperature  TQ.  The  properties  of  AT(r,z,t)  in  Eq.  (7a)  can  be  readily 

appreciated  by  discussing  two  limiting  cases  for  the  absorption 

coefficient  u-  Weak  absorption,  that  is  yL<<1 ,  means  uniform  energy 

absorption  in  z-direction.  The  z-dependent  term  in  the  second  bracket 
vanishes  and  the  remaining  factor  (1-e  uIj)  is  simply  the  fraction  of 
the  energy  absorbed  from  the  transmitted  Gaussian  beam.  Diffusion 
broadens  the  beam  profile  with  increasing  time.  Very  strong 
absorption,  that  is  yL>>1,  means  only  a  thin  layer  on  the  surface  at 
z-0  is  directly  heated  and  bulk  heat  in,  occurs  only  by  thermal 
diffusion.  This  case  is  of  interest  for  certain  practical*  dosemeter 
configurations . 
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The  actual  computer  simulations  make  use  of  the  heating  rate  of 
each  volume  element  2ir rArAz,  calculated  from  Eq.  7a: 


dT/dt 


(1-R)  a  P  e 


-r2/[4at  ♦  w2] 


.  2  I 


kLir  4at  +  w 

-an2  t 


e  m 


m-1  1  ♦  (n  /u)2 

CD 


(1-e  utjcosmir)  cos(n  z)}  (7b) 


m 


The  four  first-order  differential  equations  which  describe  the 
independent  processes  of  peaks  2  to  5  [20]  were  solved  together  with 
the  heating  rate  of  Eq.  7b  by  the  Runge-Kutta  method  for  each  cell 
Ar-Az.  Sufficient  precision  was  obtained  with  Ar  -  2w/50,  Az  »  L/50 
and  At  -  10~J  x  the  total  exposure  time.  The  sum  of  the  relative  peak 
heights  [20]  was  used  a  measure  of  the  total  trapped  electron 
population  and  each  peak  height  was  normalized  by  this  sum  to  provide 
the  initial  conditions  for  the  TL  differential  equations.  Annealing  of 
the  lowest  temperature  peak  it  2  or  any  other  is  easily  simulated  by 
putting  its  relative  height  to  zero. 


3.  EXPERIMENTAL  PROCEDURE 

The  facility  shown  schematically  in  Figure  3  was  constructed  for 
the  measurement  of  laser  stimulated  thermoluminescence  response  curves. 

A  radio-frequency  excited  C02  laser  (Laakmann  Electro-Optics, 
rated  at  4  Watts)  delivers  a  ew  beam  of  1.3mm  full  width  at 
1/e-maximum.  Together  with  all  the  required  optical  elements  and 
positioners,  this  laser  is  mounted  on  a  sturdy  aluminum  table. 

The  electromechanical  shutter  is  the  first  optical  element  in  the 
laser  beam  path.  To  avoid  thermal  distortion  of  the  shutter  blades, 
they  are  gold  plated  for  high  reflectivity.  During  the  closed  position 
of  the  3hutter,  the  laser  beam  is  reflected  from  the  blades  into  a  beam 
dump  which  simply  absorbs  and  dissipates  the  laser  power.  The  shutter 
has  a  DC  power  supply  which  is  housed  inside  the  electronic  control 
enclosure.  The  electronic  controlled  shutter  is  open  only  to  heat  a 
sample  on  the  sample  holder  or  to  measure  by  a  powermeter  the  laser 
power  delivered  to  the  sample.  The  powermeter  is  mounted  on  a  shaft 
which  slides  through  the  light  tight  box. 
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Following  the  laser  beam  path  beyond  the  shutter,  a  5x  beam 
expander  is  mounted  on  the  reader  table.  It  can  be  adjusted  to  produce 
various  beam  diameters  (1.5  to  6mm  full  width)  entering  the  focusing 
cube.  The  focusing  cube  deflects  the  beam  by  90  degrees  through  an  AR 
coated  germanium  window  and  causes  it  to  converge  on  the  TL  sample  in 
the  light  tight  box.  The  germanium  window  is  transparent  at  the  10.6pm 
laser  wavelength,  but  opaque  to  visible  light.  The  beam  position  on 
the  TL  sample  is  determined  by  moving  the  sample  holder  with  two 
connected  x-y  translation  stages.  The  translation  stages  have  two 
actuators  which  are  driven  by  the  Newport  Corporation  model  350M  motor 
drive  system  and  give  an  accuracy  of  the  beam  position  up  to  0.1pm. 

The  TL  sample  is  heated  as  it  absorbs  this  infra-red  radiation 
from  the  laser.  The  high  power  density  of  the  focused  laser  beam 
causes  the  selected  spot  on  the  TL  sample  to  be  heated  rapidly,  while 
the  reminder  of  the  sample  remains  cool. 

The  TL  emitted  by  the  sample  passes  the  sample  substrate  and 
reaches  into  the  photomultiplier  tube  (PMT)  below.  The  PMT  is 
protected  from  two  types  of  accidental  exposure,  one  due  to  the  laser 
by  a  glass  window  which  will  absorb  the  10.6pm  laser  beam,  another  one 
due  to  the  visible  light  (when  the  sample  chamber  is  open)  by  a  thick 
aluminum  slide.  The  PMT  is  an  EMI  type  9924B  with  a  1"  diameter 
bialkali  photocathode.  The  current  signal  from  the  PMT  is  digitized, 
stored  and  displayed  by  a  Data  Precision  6000  Universal  Waveform 
Analyzer  (transient  digitizer). 

The  laser  TLD  reader  was  designed  to  permit  examination  of  a 
number  of  experimental  parameters  related  to  laser  heating  TLD,  e.g., 
the  laser  power  density  that  the  sample  receives,  and  the  beam  diameter 
which  13  set  by  mechanically  adjusting  the  beam  expander.  Laser  power 
and  heating  time  are  controlled  by  settings  in  the  electronic  control 
unit  (Fig.  4). 

The  timing  circuits  in  the  electronic  control  unit  consist  of  four 
LM555  timers  that  generate  four  timing  pulses.  The  first  pulse  drives 


the  cycle  indicator  and  limits  read  cycles  to  one  per  second.  The 
second  pulse  is  used  to  turn  the  laser  off  and  open  the  shutter.  The 
third  pulse  turns  the  laser  on  for  a  precise  preselected  timing 
Interval.  The  transient  digitizer  is  triggered  at  the  beginning  of 
this  pulse,  and  the  laser  is  turned  off  at  its  end.  The  fourth  pulse 
triggers  the  closing  mechanism  of  the  shutter.  The  laser  power  during 
the  third  timing  pulse  is  determined  by  modulation  control.  Whenever 
all  four  of  these  pulses  are  off,  the  laser  operates  in  a  cw  mode  at 
maximum  output  power  to  maintain  a  stable  operating  temperature. 

The  modulation  circuit  consists  of  a  10  MHz  Oscillator  that 
triggers  a  variable  length  pulse  generator.  The  modulation  control  is 
used  to  vary  the  pulse  width  of  the  resulting  pulse  train,  permitting 
the  laser  output  power  to  be  varied  over  a  range  from  zero  to  full 
power.  The  modulation  driver  amplifies  the  modulation  signal  to  the  6 
volts  pulse  height  required  for  the  laser  RF  power  supply.  The  shutter 
control  signal  and  the  data  signal  are  both  logic  level  signals.  The 
74LS  series  TTL,  logic  circuits  are  employed  for  all  control  logic. 

As  will  be  demonstrated  below  (see  for  example  Fig.  12),  the  shape 
of  the  thermoluminescence  response  curve  depends  critically  on  the 
laser  beam  profile  in  the  sample  plane  at  z-0.  For  this  reason  we  have 
given  its  precise  measurement  particular  attention,  using  the  so-called 
scanning  knife  edge  method  [33]. 

A  photodetector  measures  the  laser  power  not  occluded  by  a  knife 
edge,  scanning  through  the  beam,  as  a  function  of  position.  The 
spatial  profile  is  then  given  by  the  first  derivative  of  the 
photodetector  signal  with  respect  to  the  scan  direction.  A  necessary 
requirement  is  that  the  mathematical  form  of  the  profile  be  separable 
in  the  coordinates  perpendicular  to,  and  parallel  to  the  scan 
directions  as  in  the  case  of  a  Gaussian. 

To  check  for  deviations  from  circular  symmetry,  we  have  always 
scanned  in  two  perpendicular  directions.  For  cw  lasers  this  procedure 
is  readily  executed  by  scanning  at  a  constant  rate  and  analog 


differentiation  with  an  appropriate  RC  circuit  [3^].  However,  for 
infra-red  lasers  such  as  the  one  used  in  our  work,  thermal  detectors 
are  commonly  employed  to  monitor  the  beam  power.  Because  of  their 
saturation  characteristics  they  are  generally  not  suitable  for  scanning 
knife-edge  applications.  Modifications  of  this  method  are 
required.  They  are  shown  schematically  in  Fig.  5.  The  cw  beam  is 
chopped  with  the  electro-mechanical  shutter  [3]  driven  by  an 
oscillator.  The  oscillator  output  is  used  simultaneously  as  a  clock 
for  a  transient  digitizer  (Data  Precision  model  DP  6000)  to  synchronize 
the  sampling  interval  with  the  occurance  of  the  signal  pulse  from  a 
piezoelectric  detector  (Barnes  Engineering  model  350  PZT).  In  our  case 
the  shutter  is  opened  every  73- 5ms  for  3.5ms.  The  recorded  detector 
signal  during  a  complete  scan  has  the  appearance  of  a  digitally 
recorded  cw  signal,  exhibiting  only  the  slow  temporal  structure  due  to 
increasing  occlusion  of  the  beam  by  the  knife  edge.  This  signal  is 
then  differentiated  mathematically  by  the  microprocessor  of  the  DP 
6000. 

Examples  of  beam  profiles  obtained  in  this  fashion  are  shown 
Fig.  6.  The  combination  of  the  beam  expander  and  5”  focal  length 
focusing  cube  afforded  variations  of  the  laser  spot  size  in  the  sample 
plane  between  0.3mm  and  2.22mm.  The  long  focal  length  of  the  cube  and 
the  maximum  laser  beam  penetration  depth  in  our  samples  of  less  than 
300u®  implies  that  the  laser  beam  profile  is  independent  of  the 
z-direction. 

In  the  following  Sections  we  present  experimental  results  obtained 
with  a  number  of  different  dosemeter  configurations  and  comparisons 
with  theoretical  calculations. 

Case  1:  LiF : Ti ;Mg  Chip 

Laser  heating  of  microcrystalline  hot-pressed  TLD-100  dosemeter 
material  produced  by  the  Harshaw/Filtrol  Corporation  [ll]  is  of 
interest  as  it  is  widely  employed  in  routine  personnel  dosimetry.  From 
a  theoretical  point  of  view  this  case  is  important  because  of  the 
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rather  modest  absorption  coefficient  (u-40  cm“*)  of  LiF  for  10.6uo 
photons.  The  laser  beam  penetrates  250un  (1/e  attenuation  depth)  into 
the  900»m  thick  chip.  As  it  turns  out,  none  of  the  previously 
reported  solutions  of  the  heat  flow  equation  (Eq.  3)  covers  this 
situation  in  conjunction  with  a  Gaussian  heating  beam. 

Unfortunately  only  small  chips  (3mm  x  3®o  x  0.9mm)  are 
commercially  available  and,  therefore,  care  has  to  be  taken  not  to 
violate  the  assmuption  of  a  semi-infinite  medium  (see  Section  II)  when 
comparing  computed  with  experimental  results.  In  practice  this  can  be 
done  by  choosing  a  beam  diameter  that  is  much  smaller  than  the  width  of 
the  sample  and  by  limiting  the  total  exposure  of  the  laser  heating  beam 
to  a  time  during  which  the  sample  edges  have  not  yet  been  heated  by 
thermal  diffusion.  In  order  to  reach,  during  this  time,  in  the  center 
of  the  exposed  area  the  highest  temperature  peak  (peak  #  5),  the  power 
density  in  the  beam  center  must  exceed  a  minimum  value.  In  essence 
this  requirement  means  that  a  substantial  part  of  the  chip  volume  must 
be  heated  to  above  the  maximum  temperature  for  which  TL  emission  will 
still  occur  before  the  edge  areas  experience  any  temperature  increase 
by  thermal  diffusion.  We  have  performed  these  experiments  successfully 
with  beam  diameters  around  0.084cm  and  approximately  5  Watt  power. 
Under  these  conditons  a  300ras  exposure  time  raised  the  temperature  at 
the  edges  by  only  a  few  degrees  K  above  room  temperature.  Continued 
exposure  to  the  .heating  beam  beyond  this  time  resulted  in  heat  pile-up 
in  the  chip  and  eventual  increase  of  the  temperature  of  the  total 
sample  to  above  700K,  a  situation  not  covered  by  the  theory  presented 
in  Section  II. 2.  However,  these  experiments  are  perhaps  important  as 
they  indicate  that  entire  TLD-100  or  similar  chips  can  indeed  be  heated 
with  an  appropriate  laser  beam  so  as  to  measure  the  total 
thermoluminescence  emission.  In  fact  we  have  been  able  to  heat  these 
chips  to  600K  within  about  2  seconds  with  a  2.5nyn  diameter  beam  of 
approximately  8  Watt  power. 

The  temperature  distribution  established  at  the  end  of  the 
exposure  (300ms)  of  a  LiF  chip  to  a  0.084cm  full  width  at  1/e-maximum 
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of  4.93  Watts  power  is  depicted  in  Fig.  7.  Isotherms,  calculated  from 
Eq.  7a  are  shown  for  temperatures  corresponding  to  the  peak 
temperatures  of  peaks  2  through  5  (383K,  421 K,  457K,  and  483K, 
respectively).  The  volume  inside  the  483K- isotherm  is  heated  to 
temperatures  for  which  thermoluminescence  emission  from  peaks  2  through 
4  has  ceased  at  300ms.  Fig.  7  illustrates  the  nonuniformity  of  both 
the  temperature  distribution  and,  albeit  indirectly,  the 
thermoluminescence  emission  pattern. 

We  have  made  no  attempts  in  our  experiments  to  measure  the  spatial 
distribution  of  the  time-dependent  thermoluminescence  response  of  the 
TLD-100  chip  to  the  Gaussian  heating  beam.  Instead  the  total  emission 
IfL ( t )  was  monitored  with  a  photomultiplier  tube  and  compared  with 
computations  performed  with  the  theory  outlined  in  Section  II. 2.  The 
results  are  presented  in  Fig.  8  for  a  TLD-100  chip  that  was 
pre-annealed  to  373K  for  ten  minutes  in  order  to  remove  peak  2.  This 
type  of  thermal  treatment  after  exposure  to  the  ionizing  radiation  and 
before  measuring  the  thermoluminescence  glow  curve  is  customary  in 
dosimetry  applications  of  this  material.  The  agreement  between  the 
calculated  and  experimental  glow  curve  is  rather  satisfying  in  light  of 
the  uncertainties  concerning  the  exact  electron-kinetic  model 
description  of  the  thermoluminescence  emission  from  this  material  and 
the  fact  that  the  TLD-100  chip  is  not  a  single  crystal  as  assumed  in 
the  calculations.  An  additional  uncertainty  stems  from  the  rather 
strong  temperature  dependence  of  the  thermal  conductivity,  k. 
According  to  Men  et  al .  [35]  k  decreases  with  increasing  temperature 
approximately  as  T“ 1 .  The  calculated  curve  shown  in  Fig.  8  was 
obtained  with  k-0.0275  W^/cmK,  which  is  approximately  one  third  of  the 
value  measured  for  a  single  crystal  at  room  temperature.  We  have  also 
performed  the  calculation  with  an  exact  solution  of  the  heat  flow 
equation  with  k(T)-kQ  Tq/T  using  Kirchhoff  transformation  techniques 
(see  Appendix  A).  However,  no  significant  improvement  in  the  agreement 
between  measured  and  calculated  thermoluminescence  response  curves  was 
obtained,  perhaps  indicating  that  unaccounted  for  experimental  factors 
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lnfluenced  the  results  more  than  an  improved  theory  of  heat  conduction. 
We  suspect  that  the  mi  cocrystalline  "ceramic"  nature  of  the  sample  and 
the  roughness  of  the  unpolished  surface  caused  some  deviations  from  the 
implicit  assumption  that  the  optical  and  thermal  properties  of  this 
material  are  very  close  to  those  of  single  crystalline  LiF.  Scattering 
of  the  laser  light  probably  results  in  some  broadening  of  the  Gaussian 
beam  profile  Inside  the  chip,  and  the  thermal  conductivity  may  be 
lower  than  that  of  a  single  crystal.  In  addition,  it  is  well  known 
that,  in  order  to  reproduce  the  exact  relative  peak  heights  of  peaks  2 
to  5  of  the  TLD-100  phosphor  as  used  in  ref.  [20],  the  annealing 
procedure  must  be  reproduced  exactly  [ 1 2 ] .  Slight  deviations  might 
have  been  present  as  compared  to  those  used  by  McKeever  [ 20 ]  ,  resulting 
in  small  peak  height  differences  particularly  of  the  low  temperature 
peaks  which  significantly  influence  the  thermoluminescence  response 
curve  obtained  by  Gaussian  beam  heating  (see  Case  4  of  this  Section). 


Case  2:  Polyimide  foil  Loaded  with  LiF : Mg ,Ti  Powder 

Thin  dosemeters  are  desirable  for  applications  in  the  dosimetry  of 
low  energy  beta  rays.  We  have  experimented  with  84ym  thick  foils  of 
polyimide  (Kapton,  DuPont)  which  were  produced  by  thoroughly  mixing 
TLD-100  powder  (grain  size  20  to  30yffl)  with  liquid  polyimide  and 
spreading  this  slurry  uniformly  onto  glass  slides.  After  curing  for 
several  minutes  at  373K  and  subsequently  for  1  hour  at  573K,  the 
translucent  yellow  film  was  removed  by  immersion  in  water.  Correcting 
for  the  weight  loss  during  the  polymerization  process,  the 
powder- polyimide  weight  ratio  of  the  dry  foil  was  about  2:1. 

Pieces  of  2.5  x  2.5cm2  size  were  mounted  in  metal  frames  and 
exposed  to  Y-rays  from  a  *°Co  source.  The  absorbed  dose  was  26  Gy. 
Following  a  10  minute  pre-anneal  at  373K,  the  thermoluminescence 
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response  curves  were  measured  for  a  fixed  laser  power  of  4.78  Watts  and 

0.0^2 

various  laser  spot  sizes  ranging  from  0.039  to  0>>r€cra  diameters.  A 
family  of  curves  measured  under  these  conditions  is  depicted  in  Fig.  9* 
Several  different  sites  on  the  same  sample  were  heated  for  each  laser 
spot  size  to  assess  sample  uniformity  and  laser  stability. 

All  experiments  were  performed  to  test  the  validity  of  the  theory 
presented  in  Section  II. 2  rather  than  to  fabricate  viable  dosemeters. 
For  example,  Kapton  absorbs  quite  strongly  the  TLD-100 
thermoluminescence  emission  centered  around  400nm  and,  therefore,  a 
foil  of  this  type  is  rather  insensitive  as  a  dosemeter.  Nevertheless, 
important  information  on  the  energy  transfer  from  a  laser  beam  to  a 
sample  of  this  type  and  the  resulting  thermoluminescence  response  can 
be  gathered  and  eventually  applied  to  Kapton  foils  loaded  with  TL 
powders  other  than  LiF. 

The  TLD-100  loaded  Kapton  foil  does  not  strongly  absorb  10.6um 
photons.  Yet,  direct  measurement  of  its  absorption  coefficient  is  not 
trivial  because  of  the  photon  scattering  off  the  LiF  particles 
dispersed  in  it.  Since  the  foil  is  quite  thin  and  LiF  ha3  a  rather 
small  absorption  coefficient  for  10.6um  photons,  one  can  safely  assume 
that  the  heat  deposition  is  independent  of  z  throughout  the  entire 
thickness  L.  However,  scattering  introduces  an  element  of  uncertainty 
whose  influence  on  the  agreement  between  theory  and  experiment  can  only 
be  estimated  by  the  comparison  of  measured  and  calculated 
thermoluminescence  response  curves.  This  case  of  relatively  weak 
absorption  is  characterized  by  pL<<1.  As  a  result,  Eq.  7a  becomes 
independent  of  z: 


t  -r2/[4a(t-x)  +  w*] 

AT  -  B  /  dx  ® -  (8) 

t-o  4a(t-t)  +  w1 


with 


-20- 
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(9) 


In  attempting  to  calculate  the  family  of  curves  in  Fig.  9a,  the 
problem  of  determining  the  thermal  and  optical  properties  of  the  sample 
that  enter  Eq.  8  presents  itself.  For  example,  it  is  not  immediately 
possible  to  obtain  the  relevant  averaged  thermal  conductivity,  k,  for 
such  a  heterogeneous  mixture  from  the  known  values  of  its  components. 
Also,  the  total  absorbed  power  P(l-R)(1-e  w1j)  is  difficult  to  determine 
without  a  suitable  calorimeter,  while  the  averaged  specific  density,  p, 
and  heat  capacity,  c,  can  be  calculated  in  principle  from  the  values  of 
the  components. 

However  the  quantity  B  defined  by  the  entire  parameter  combination 
in  Eq.  (9)  can  be  measured  by  exploiting  the  normally  undesired  onset 
of  incandescent  emission  that  is  observed  when  the  sample  is  exposed  to 
the  laser  beam  significantly  beyond  the  occurrence  of  the  glow  curve 
maximum  (see  Fig.  9b).  We  have  developed  a  method  for  this  task  that 
is  based  on  Planck's  law  of  blackbody  emission.  It  makes  use  of  the 
exponential  time  dependence  of  thi3  emission  when  it  Is  measured  for  a 
narrow  wavelength  Interval  whose  center  wavelength,  X,  has  to  be 
selected  in  the  red  or  better  yet  in  the  infra-red  so  as  to  not  destroy 
the  sample  by  overheating  it.  Also,  the  experiment  has  to  be  performed 
with  significantly  higher  laser  power  densities  as  compared  to  the  ones 
used  in  the  thermoluminescence  experiments  in  order  to  assure  that  the 
condition  Hat<<w2  is  fulfilled  so  that  measurable  incandescent  emission 
is  generated  before  the  onset  of  thermal  diffusion.  Under  this 
condition  the  temperature  rise  of  the  heated  spot  is  linear  in  time  as 
can  be  seen  from  Eq.  8. 
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and  It  can  be  shown  (see  Appendix  B)  that  the  so  measured  incandescent 
emission  intensity,  I,  obeys  the  relation 

(|£)/I  -  [a/ (Bt/w*  ♦  T0)a](B/ws)  (11) 

where  a  -  hc/Xk  and  TQ  is  the  Initial  temperature.  The  quantity  B  earn 
be  calculated  by  carefully  measuring  the  laser  spot  size,  2w,  and  dl/I 
as  a  function  of  time.  Details  of  this  method  are  presented  in 
Appendix  B. 

For  the  above  mentioned  TLD-1 00-loaded  Kapton  foils  B/P  was  found 
to  be  5-53  [Kcm2/j].  With  this  knowledge  and  the  known  laser  power  the 
thermoluminescence  response  curves  of  Fig.  9a  can  be  calculated  from 
Eq.  8  and  the  computational  procedure  for  glow  curves  obtained  with 
locally  non-uniform  heating  as  described  in  Section  II.  This  was  done 
by  normalizing  a  single  member  of  a  set  of  curves  measured  for 
different  laser  beam  diameters,  2w,  and  selecting  a  value  for  the 
unknown  thermal  diffusivity,  a,  30  as  to  obtain  the  best  fit  for  all 
members  of  the  set.  Usually  a  rather  good  fit  is  possible  for  the 
curves  measured  with  large  power  densities  of  the  laser  beam  (e.g.  in 
the  present  case  smaller  spot  sizes  for  the  fixed  laser  power  of  4.70 
Watts).  Larger  spot  sizes  tend  to  yield  computed  curves  that  rise 
later  than  the  experimental  ones  and  are  steeper  prior  to  reaching  the 
maximum.  Perhaps  the  scattering  of  the  laser  light  from  the  LiF 
grains,  disregarded  in  our  computations,  and  the  associated  distortion 
of  the  laser  beam  profile  may  be  a  reason  for  this  discrepancy.  The 
value  obtained  for  the  thermal  diffusivity  a-0.0175  cmVs  appears 
reasonable  In  light  of  the  values  for  LiF  (0.02743cm2/s)  and  Kapton 
(O.OOlcmVs). 

In  closing  this  discussion  of  the  powder-loaded  Kapton  foils  it 
may  be  of  interest  to  mention  that  the  measurement  of  laser-heated  TL 
glow  curves  and  the  developed  computational  procedure  presented  in  this 
work 


l 
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polnt  to  a  new  method  for  the  determination  of  the  thermal  properties 
of  otherwise  difficult  to  measure  thin-foil  TLD  dosemeters. 

Having  completely  determined  the  thermal  and  optical  properties  of 

a  given  dosemeter  of  this  type,  it  is  relatively  straightforward  to 

simulate  its  thermoluoinescence  response  as  a  function  of  other 

experimental  parameters  such  as  laser  power,  beam  diameter,  foil 

thickness,  pre-annealing  temperature  and  duration  and  the  like.  As  an 

example  Fig.  10  shows  the  dependence  on  the  laser  power  of  the  sample 

0>  0 

from  Fig.  9a  for  a  fixed  beam  diameter  2w-0^£fcm.  Thus,  the  techniques 
are  now  available-  for  the  optimization  of  dosemeter  designs  for  the  use 
in  laser-heated  thermoluminescence  dosimetry  readers  [3]. 

Case  3:  Thin  LiF:Mg,Ti  Phosphor  Layer  on  a  Glass  Slab  -  a)  Heating 
Beam  on  Layer 

Another  dosemeter  configuration  of  potential  importance  for 
laser-heated  thermoluminescence  dosimetry  is  a  thin  TL  phosphor  layer 
on  a  substrate  that  is  transparent  for  the  emitted  thermoluminescence 
photons.  In  this  arrangement  the  heating  beam  may  be  directed  so  as  to 
expose  a  small  spot  of  the  sample  on  the  side  that  is  covered  with  the 
active  layer  and  the  thermoluminescence  emission  is  monitored  with  a 
photon  detector  on  the  opposite  3ide.  The  substrate  may  be  chosen 
according  to  its  absorption  coefficient  for  the  photons  of  the  heating 
beam.  For  example  BaF,  LiF,  and  sapphire  or  glass  are  materials  having 
absorption  coefficients  increasing  in  that  order.  Assuming  the 
phosphor  layer  is  sufficiently  thin  and  composed  mainly  of 
TLD- 100  phosphor  powder,  little  energy  is  transferred  from  the  10.6vim 
beam  U3ing  a  BaF  substrate.  However,  useful  thermoluminescence 
response  curves  are  obtained  with  a  single  crystalline  LiF  substrate 
(see  Fig.  11).  The  most  efficient  utilization  of  the  laser  beam  energy 
is  achieved  by  a  highly  absorbent  substrate  such  as  common  boroslicate 
glass.  The  thermoluminescence  response  of  3uch  a  dosemeter 
conf iguration  is  the  rather  complex  result  of  both  direct  and  indirect 
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(via  absorption  by  the  glass  and  subsequent  diffusion)  heating  of  the 
phosphor.  Both,  its  optical  properties  at  10.6um  a3  well  as  its 
thermal  characterstics,  are  due  to  the  contributions  of  the  individual 
properties  of  both  components  In  an  unknown  fashion  and  are  for  this 
reason  difficult  to  measure.  However,  by  applying  a  layer  that  is 
considerably  thinner  than  the  substrate  and  considering  the  fact  that 
this  layer  does  only  moderately  absorb  the  laser  photons,  heating  is 
mostly  indirect  via  diffusion  from  the  strongly  absorbing  glass 
surface.  Therefore,  one  can  assume  that  the  thermal  diffusivity,  a,  of 
the  dosemeter  is  very  close  to  that  of  the  glass  (a  =*  0.0075  cm2/s). 
In  fact,  it  should  be  slightly  higher  than  that  of  lime  glass  because 
of  the  contribution  to  the  "conf igurat ion-averaged"  a  from  the  LiF 
content  of  the  phosphor  layer.  We  have  performed  all  calculations  (see 
below)  of  the  thermoluminescence  reponses  of  thin-layer  dosemeter 
configurations  which  use  borosilicate  glass  substrates  with  the  value 
of  a  -  0.01  cm2/s  (glass  composition:  64t  silicon  dioxide,  1 3%  alkali 
oxides,  8t  boron  oxide,  and  1 5%  other  oxides).  The  only  parameter 
combination  not  known  for  this  type  of  dosemeters  is  B  (see  Eq.  9), 
which  is  obtained  by  the  best  fit  of  calculated  and  measured  TL 
response  curves. 

We  have  performed  experiments  with  thin  uniform  TLD-’OO  powder 
layers  in  a  silicone  binder  on  0.1 59mm  thick  borosilicate  glass 
microscope  cover  slides.  An  example  of  such  a  family  of  glow  curves 
measured  with  a  fixed  laser  power  and  for  various  laser  spot  sizes  is 
presented  in  Fig.  12.  A  25  x  25mm2  slide  was  coated  with  a  51  ym  thick 
layer  of  LiF  (TLD-100)  powder  (grain  size  between  30  and  40ym)  in  a 
silicone  binder  (cured  weight  ratio  2:1)  and  exposed  to  50  keV  x-rays. 
After  10  rain,  pre-anneal  at  373K  to  remove  peak  2  [ 20 ] ,  well  separated 
areas  on  thi3  sample  were  exposed  to  laser  heating  pulses  of  4.45  Watt 
power  and  spot  sizes  of  full  widths  ranging  from  1  to  2mm.  In  order  to 
avoid  the  onset  of  Incandescent  emission  from  the  center  of  the  exposed 
spot  the  exposure  times  were  adjusted  according  to  the  power  density, 
ranging  from  100ms  to  250ms.  From  the  best  fit  of  the  computed  curves 
to  these  measurements  B  was  found  to  be  20  KcraVs  for  the  laser  power 
stated  above. 


With  this  value  of  B  other  aeries  of  thermoluminescence  response 
curves  for  this  type  of  dosemeter  can  readily  be  generated.  For 
example  the  influence  of  the  laser  power  P  for  various  spot  sizes  or 
the  thickness  L  of  the  sample  can  be  studied  by  computer  simulation  in 
order  to  arrive  at  an  optimal  dosemeter  design  for  a  given  available 
COj  laser.  An  example  is  presented  in  Fig.  12b  which  demonstrates 
rather  large  changes  in  the  peak  location  of  the  thermoluminescence 
response  curves  when  varying  the  laser  power  by  only  10t,  indicating 
the  necessity  to  use  a  highly  stabilised  laser  beam  for  applications  of 
laser  heating  in  thermoluminescene  dosimetry. 

Case  4:  Thin  LiF:Mg,Tl  Phosphor  Layer  on  a  Glass  Slab  -  b)  Heating  Beam 
on  Glass  Side  Opposite  the  Layer 

It  is  of  course  also  possible  to  heat  the  active  thermolumiscent 

layer  through  the  substrate.  Choosing  again  a  relatively  thin 

microscope  cover  slide  interesting  thermolumine3cence  response  curves 

from  a  TLD-100  phosphor  layer  were  obtained.  Compared  to  those 

measured  by  frontal  heating  with  the  same  laser  power  and  beam 

diameter,  they  all  peak  at  a  later  time  because  diffusion  through  the 

strongly  absorbing  substrate  in  now  the  only  mechanism  of  heat 

transfer.  At  the  same  time  the  effected  spot  size  is  considerably 

larger.  ^ 

A  series  of  curves  measured  with  a  laser  power  of  JJtdT  Watts  and 

1.6T  a.|  3 

spot  sizes  ranging  fr'om  to  mm  full  width  is  shown  in  Fig.  13 

for  the  same  sample  used  in  the  experiments  of  Fig.  12.  The 

theoretical  curves  were  again  computed  from  Eq.  7b  with  an 

experimentally  determined  value  for  the  parameter  combination 

B-16.75  cm2/s.  The  best  fit  was  obtained  with  a  thermal  diffusivity  of 

a-(7.8  t  0.5)  x  10_2cm2/s.  The  calculated  dependence  of  the 

i.f>r 

thermoluminescence  response  on  the  laser  power  for  a  spot  size  of  J^Stnm 
Is  plotted  in  Fig.  14  together  with  curves  calculated  for  a  fixed  power 
and  3 pot  size  and  different  sample  thickness  L. 
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We  have  used  the  same  sample  in  an  attempt  to  study  the  effect  of 
pre-annealing  on  the  thermoluminescence  response  curve.  These 
experiments  are  not  just  a  trivial  repetition  of  what  has  been  done 
before  by  many  authors  to  find  the  appropriate  treatment  for  removal  of 
rapidly  fading  low  temperature  glow  peaks  prior  to  the  dose 
measurement.  Because  the  temperature  distribution  is  nonuniform  when  a 
semi-infinite  slab  is  heated  by  a  localized  beam  of  Gaussian  profile, 
uniform  pre-annealing  effects  the  response  curve  much  more  severely 
than  the  conventional  glow  curve.  This  becomes  obvious  by  careful 
examination  of  the  photographs  of  Fig.  2  which  may  be  considered  as  a 
"coded"  temperature  distribution  in  the  emitting  layer  for  the  various 
times  after  onset  of  the  laser  exposure.  Imagine  ring-like  surface 
areas  between  two  adjacent  isotherms  of  a  given  small  temperature 
difference.  At  the  perimeter  of  the  emission  pattern,  which  is  at  a 
low  temperature,  the  emission  from  such  a  ring  is  due  to  the  low 
temperature  peaks  of  the  conventional  glow  curve.  A  ring  of  the  same 
temperature  difference,  but  located  closer  to  the  center,  emits  at  the 
same  time  at  a  higher  temperature.  Assuming  for  the  3ake  of  simplicity 
identical  brightness  of  these  two  areas,  their  contribution  to  the 
total  emission  (monitored  with  a  detector  and  without  spatial 
resolution)  is  proportional  to  the  square  of  their  respective  radii. 
Removing  the  low  temperature  glowpeaks  by  uniform  pre-anneal  treatment 
affects  the  outer  and,  thus,  larger  emitting  area  and  consequently  the 
strongest  contributions  to  the  emission  at  that  time.  Fig.  15  is  an 
illustration  of  this  phenomenon.  Three  thermoluminescence  response 
curves  are  3hown  measured  with,  a  laser  power  of  3.63  Watts  and  a  spot 
size  of  1.6mm  1 /e-power  diameter.  The  sample  is  the  3ame  as  that  used 
to  produce  the  curves  of  Fig.  13.  A  373K  pre-anneal  of  only  5  min., 
which  is  known  to  remove  the  low  intensity  peak  2  of  the  conventional 
glow  curve  [ 20  ]  has  a  rather  large  effect  on  the  thermoluminescence 
response  curve  measured  with  Gaussian  laser  beam  heating. 
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Pre-anneallng  at  the  same  temperature  for  10  min.  has  almost  no 
measurable  effect  on  peaks,  3.  ^  and  5  of  the  conventional  glow  curve 
while  the  difference  can  still  be  detected  in  curve  c  of  Fig.  15.  We 
refrain  here  from  reporting  pre-annealing  measurements  and  computer 
simulations  performed  with  all  the  sample  types  used  in  our 
investigations.  Instead  we  restrict  ourselves  to  pointing  out  that  in 
practical  dosimetry  applications  of  laser  heating  continuous  dosemeter 
layers  with  Gaussian  beam  profiles  [15]  great  care  has  to  be  taken  to 
reproducibly  pre-anneal  the  samples.  In  addition,  fading  of  the 
Information  on  absorbed  loses  stored  in  the  thermoluminescence 
dosemeter  will  cause  a  similarly  enhanced  error  as  compared  to 
conventional  contact  heating. 

III.  Conclusions 

We  have  presented  the  theory  of  thermoluminescence  response  curves 
generated  with  Gaussian  laser  beam  profiles  for  a  number  of  practical 
and  interesting  dosemeter  configurations,  and  we  have  verified  our 
theoretical  approach  by  direct  comparison  with  experiment.  We  have 
also  shown  that  by  computer  simulation  one  can  assess  the  effects  of 
such  design  parameters  as  layer  thickness,  laser  power  and  beam  size, 
pre-annealing  temperature  and  duration,  substrate  material  and  its 
thickness. 

The  extremely  fast  heating  rates  (>  iC'K/s)  which  can  be  obtained 
by  laser  heating  versus  that  of  conventional  contact  heating  in 
solid-state  dosimetry  (-  IK/s)  and  the  reduced  amount  of  phosphor 
required  in  our  thin  film  configuration  indicate  quite  clearly  that  the 
techniques  which  we  have  developed  and  presented  here- in  while  of 
interest  in  their  own  right  will  be  of  significant  practical 
importance. 
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Figure  Captions 

Fig.  1  Schematic  representation  of  a  Gaussian  laser  beam  profile 
exposing  the  surface  (z  -  0)  of  a  slab  of  thickness  L. 

Fig.  2  Thermoluminescence  emission  pattern  obtained  from  a  50  yin  thick 
uniform  ZnS:Cu  layer  on  a  1  mm  thick  glass  slab  (25  x  25  mm2).  The 
sample  was  heated  by  a  cw  C0a  laser  beam  of  4.3  Watts  and  5.4  mm  full 
width  at  1/e  peak  power.  The  beam  was  directed  perpendicularly  at  the 
center  of  the  glass  slab  on  the  side  opposite  to  the  phosphor  layer. 
The  photographs  were  taken  at  the  indicated  times  after  onset  of  the 
laser  exposure. 

Fig.  3  Schematic  representation  of  the  laser-heated  thermoluminescence 
reader . 

Fig.  4  Block  diagram  of  the  electronic  control  system  for  the 
apparatus  shown  in  Fig.  3. 

Fig.  5  Schematic  arrangement  for  knife-edge  scanning  of  the  C02  laser 
beam  profile.  Mote  that  the  position  of  the  knife-edge  in  the  sample 
plane  is  shown  here  vertically  because  the  focusing  cube  (see  Fig.  3) 
is  omitted  for  simplicity.- 

Fig.  6  Example  of  a  typical  nearly  Gaussian  spatial  profile  obtained 
for  the  CO*  laser  beam  at  z  -  0  (see  Fig.  1). 

Fig.  7  Temperature  distribution  established  after  heating  a  3mall  LiF 
slab  of  area  3x3nim2  and  thickness  L-0.9mm  for  300ms.  A  Gaussian  laser 
beam  of  4.93  Watts  power  and  .84ram  full  width  impinges  on  the. center  of 
the  slab  at  z/L-0.  The  calculations  were  performed  with  Eq.  7b  and 
k«0.0275  W/cmK,  c-1.56  Ws/gK,  p-2.635  g/cm*,  and  an  absorption 
coefficient  for  lOum  of  y-40cm"‘. 
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Fig.  8^  Measured  ^and  calculated  (solid  line)  thermoluminescene  response 
curves  obtained  from  the  LiF:Mg,  Ti  (TLD-100)  chip  of  Fig.  7  after 


x-ray  exposure.  The  calculations  were  performed  with  Eq.  7b  and  the 


first-order  electron  kinetic  parameters  (frequency  factor,  trap  depth) 
determined  by  McKeever  [ 20 ]  for  peaks  3~5  of  this  thermoluminescent 


phosphor.  Peak  2  was  removed  by  10  min.  pre-annealing  at  373K. 


Fig.  9  a)  Measured  and  calculated  (3olid  line)  thermoluminescence 
response  curves  obtained  after  x-ray  exposure  of  a  84um  thick  polyimide 
foil  loaded  with  LiF:Mg,Ti  (TLD-100)  particles  of  20~30um  grain  size  at 
a  polymer- powder  weight  ratio  of  1:2.  The  curves  were  measured  with  a 
4.7  Watt  C02  laser  beam  of  full  widths  (from  left  to  right)  of  0.39mm, 
0.62mm,  0.77mm,  and  0.985  mm.  The  calculations  were  performed  with  Eq. 
8  and  B  -  26  Kcm2/s  (see  Eq.  9)  which  was  determined  by  the 
incandescence  emission  method  (Appendix  B)  at  this  laser  power.  The 
thermal  diffusivity  of  the  polyimide-LiF  composite  was  found  to  be 
a  -  0.0175  cm2/s  by  the  best  fit  of  the  calculation  to  a  series  of 
measured  curves  obtained  with  various  laser  spot  sizes  but  otherwise 
identical  conditions. 


Fig.  9b)  A  set  of  thermoluminescence  response  curves  obtained  with  an 
identical  sample  as  the  one  used  in  Fig.  9a  except  for  longer  heating 
times,  which  were  increased  to  demonstrate  the  occurrence  of 
incandescence  emission.  Each  curve  drops  off  sharply  upon  turning  off 
the  laser  beam. 

Laser  power:  4.75  Watts 

Laser  beam  diameter  2w  [mm]  (from  left  to  right):  0.899,  0.9*12,  0.985, 
1 .024,  1 .063,  and  1 .102. 
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Fig.  10  Example  of  calculated  dosemeter  performance.  The  sample  is 
the  self-supported  polyimide- powder  foil  from  Fig.  10,  all  optical  and 
thermal  properties  being  identical.  The  calculations  were  carried  out 
for  a  fixed  laser  spot  size  at  0.62mm  and  laser  powers  P  of  10,  5,  4,  3 
and  2  Watts.  Mote  that  for  clarity  these  thermoluminescence  response 
curves  are  scaled  by  dividing  the  TL  intensity  by  a  numerical  factor 
equal  to  the  laser  power  in  Watts. 

Fig.  11.  Thermoluminescence  response  curves  measured  from  a  uniform 
LiF :Mg,Ti  layer  (thickness:  3.7  mg/cra*,  20um  grain  3ize,  Dow  Corning 
805  silicone  binder)  on  a  0.15cm  thick  non-thermoluminescent  LiF  single 
crystal  of  25  x  25mm*  area. 

Laser  power:  P  -  4.95  Watts 
Laser  spot  size  2w(in  mm) 

1:  0.326,  2:  0.327,  3=  0.330,  4:  0.336,  5:  0.340mm. 

To  check  the  reproducibility  of  the  measurements,  several 
thermoluminescence  responses  from  different  sites  on  the  sample  were 
measured  for  otherwise  identical  conditions.  Repeated  measurements  are 
shown  here  only  for  the  largest  three  of  the  five  different  spot  sizes 
chosen  for  this  experiment. 

Fig.  12a  Thermoluminescence  response  curves  measured  with  a  4.45  Watts 
laser  beam  by  heating  the  phosphor-coated  3ide  of  a  0.159  mm  thick 
glass  slab  after  x-ray  exposure.  The  51  vim  thick  coating  consisted  of 
TLD-100  powder  (grain  size  30-40uffl)  in  a  Dow  Corning  805  silicone 
binder  (2:1  weight  ratio).  The  solid  lines  are  curves  computed  from 
Eq.  7b  and  the  first-order  electron  kinetic  parameters  determined  by 
McKeever  [ 20 ]  for  peaks  3”5.  Peak  2  was  removed  by  a  10  min. 
pre- anneal  at  373K. 
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Wlth  the  known  thermal  dlffusivity  a-0.01cm2/3  the  parameter  B  was 
determined  by  fitting  a  family  of  curves  obtained  for  the  following 
full  widths,  2w,  [mm]  (curves  from  left  to  right):  2.03,  1  .83.  1-65, 
1.45,  1.24  respectively. 

Fig.  12b  Computed  thermoluminescence  response  curves  for  various  laser 
powers.  The  sample  is  identical  to  that  of  Fig.  12a.  The  laser  spot 
size,  2w,  is  l.944mm  and  the  laser  powers  (curves  from  left  to  right) 
are  4.895,  4.45,  and  4.005  Watts,  respectively. 

Fig.  13  Thermoluminescence  response  cuves  of  the  sample  from  Fig.  12, 
however,  the  laser  heating  beam  impinges  onto  the  side  opposite  to  the 
phosphor-coated  one.  Thermal  diffusion  through  the  thin  glass  3lab  i3 
now  required  to  produce  the  thermoluminescence  emission.  With  the 
known  thermal  diffusivity  of  this  sample  (see  Fig.  12a)  the  B  parameter 
(Eq.  9)  was  determined  to  be  16.5  Kcm2/s  at  P-4.18  Watts  by  the  beat 
fit  of  computed  curves  to  the  entire  series  of  curves  measured  for 
different  laser  beam  widths  2w.  Again,  peak  2  of  TLD-100  was  removed 
by  annealing.  The  solid  lines  are  the  curves  computed  with  Eq.  7b  and 
on  the  basis  of  first-order  electron  kinetics  and  the  following 
measured  beam  widths,  2w  [mm]:  1.65,  1.994,  2.13  respectively. 


I 


! 
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Fig.  14  Thermoluminescence  response  curves  of  the  sample  from  Fig.  13 
computed  from  Eq.  7b  for  various  laser  powers  P  and  total  sample 
thicknesses  L  and  fixed  spot  size  2w-i.65mm.  Again  it  was  assumed  that 
peak  2  of  the  TLD-100  phosphor  was  removed  by  pre-annealing.  The  solid 
curves  were  calculated  for  L-0.210mm  and  the  following  laser  powers 
(from  left  to  right):  6.27,  4.18,  and  3.125  Watts. 

The  dashed  curves  demonstrate  the  effect  of  the  total  thickness, 
L,  by  increasing  the  substrate  thickness  only  (fixed  phosphor  layer 
thickness  of  0.051am).  The  calculations  were  performed  with  a  laser 
power  P—4 .18  Watts  and  the  following  total  thicknesses  L  (from  left  to 
right):  0.105,  0.210,  0.420,  and  0.630mm. 

Fig.  15  The  effect  of  pre-annealing  for  10  min.  at  373K  on  the 
thermoluminescence  response  curves  obtained  from  the  sample  of  Fig.  12 
by  heating  it  through  the  glass  with  P-3.63  Watts  and  a  laser  spot  size 
2w-1.6mm.  The  first  curve  (on  the  left)  is  obtained  without  annealing 
immediately  after  exposure  of  the  sample  to  50  kev  x-ray  for  4  min. 
The  other  two  curves  represent  the  TL  response  after  5  min.  and  10  min. 
pre-annealing  at  373K,  respectively. 

Fig.  16  Determination  of  B-(1-R)(1-e  y^)P/irpcL  (see  Eq.  9)  for 
powder-loaded  Kapton  foils  by  the  Incandescence  method  described  in 
Appendix  B.  Shown  are  the  rise  of  the  incandescence  emission,  I, 
(solid  line)  and  its  time-der  ivai  ve ,  ,31. /T,  (dashed  line)  as  a 

function  of  time  after  exposure  of  the  foil  (without  previous  exposure 

to  ionizing  radiation  to  eliminate  thermoluminescence  emission)  to  a 

4.8  Watt  laser  beam  of  3pot  size  2w-0.1l8cm  .31..  From  the  ratio 

(7t) 

,  3IW_,  determined  between  zero  and  .03s,  B  was  found  to  be  27±3  KcmVs 
7t  71 

for  this  laser  power. 


APPENDIX  A 


When  the  thermal  conductivity  varies  with  temperature,  the 
non-linear  differential  equation  of  heat  conduction  (compare  with 
Eq.  3)  can  be  written  as 

7.[k(T)VT]  ♦  g(r ,z,t )  -  pc3T/3t  ( 1 -A ) 


where  the  thermal  conductivity,  k,  and  the  thermal  diffusivity,  a,  are 
assumed  to  be  temperature  dependent;  but  the  heat  source  function, 
g(r,z,t),  does  not  depend  on  temperature. 

Equation  (1.A)  can  be  linearized  by  a  transformation  which 
involves  a  change  of  the  dependent  variable.  Ihi3  method  is  called  the 
Kirchhoff  transf orraation  [3?].  The  new  variable  of  transformation,  U, 
is  defined  as 


U 


T 

/ 

T 


k(T'  ) 


dT' 


or  by 


3U 

3T 


k(T)/ko, 


(2-A) 


where  T  is  a  reference  temperature  and  k  is  the  value  of  k(T)  at  T  . 

0  o 

Upon  substitution,  Eq.  1-A  becomes 

72U(r ,z,t )  ♦  1/kQ  g(r ,z,t )  »  1 /a  U(r ,z,t )  ,  (3"A) 

where  a  -  k(T)/pc. 

For  the  materials  of  interest  in  thermoluminescence  dosimetry,  the 
thermal  diffusivity,  <1  m  a(T),  varies  very  little  with  temperature  and 
can  be  assumed  constant.  The  initial  and  boundary  conditions  for 


« 


« 


« 


Eq.  3-A  can  be  derived  from  the  original  conditions  [9]  using  2-A,  and 
the  solution  has  the  same  form  Eq.  7a. 

For  the  case  of  LiF  [35],  the  thermal  conductivity  may  be  fitted 
by  a  function  of  the  form 


k(T)  -  Tq/T  kQ  ( 4-A) 

In  our  case,  Tq  Is  the  room  temperature.  Then  the  Kirchhoff 
transformation  gives 


U  -  Tq  Ln  (T/TQ)  (5-A) 

and  the  temperature  expression  in  terras  of  U  can  be  extracted  from 
Eq.  5~A  as 


T(r ,z,t )  -  Tq  exp(U(r,z,t)/TQ) 


(6-A) 


where  U(r,z,t)  is  given  as: 


_ _ _  U-R)oP  '  8-rVt,“<t-T)  *  “'1 

U(r,z,tj  -  — rT77 —  /  dt  -  x 


nkL 


t*o  4a(t-x)  ♦  w* 


(7-A) 


•  -am*ir*(t-t)/Ll 

2  l  - - -  (1-cos  ran )  cos(^l£)} 

m-1  1  *  L 

uL 
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APPENDIX  B 


The  parameter  combination  defining  B  in  Eq.  9  of  Section  III. 2 
contains  the  thermal  properties,  the  reflectivity,  and  the  absorption 
coefficient  for  the  TL  sample  as  well  as  the  power  of  the  laser  heating 
beam.  In  principle,  these  may  be  measured  separately.  However,  that 
turns  out  to  be  rather  difficult  for  the  self-supporting 
thermoluminescent  foils,  because  the  material  is  a  composite  of  a 
powder  and  a  polymer.  Therefore,  we  developed  a  method  to  determine  B 
experimentally  so  as  to  be  able  to  calculate  the  temperature 
distribution  T(r,z,t)  and  the  resulting  glow  curves  of  these  dosemeter 
configurations,  whose  total  absorption  of  the  laser  heating  beam  is  not 
too  large  (no  dependence  on  z).  It  will  be  shown  that  it  is  possible 
to  extract  B  from  the  incandescent  background  emission  of  the  heated 
spot  on  the  sample. 

As  described  in  Section  III,  incandescent  emission  from  the  center 
hot  spot  of  a  sample  exposed  to  the  Gaussian  laser  beam  is  always 
present  when  TL  glow  curves  are  measured  and  the  laser  beam  is  not 
prematurely  turned  off  (see  Fig.  9b). 

The  energy  density  for  a  given  wavelength  X  of  the  incandescent 
blackbody  emission  is  given  by  Planck’s  law: 

-  8TrhcX_5[exp(-^_)-l  ]~*  (1-B) 

where  T  is  the  absolute  temperature,  and  h  and  k  have  the  usual 

meaning.  We  tried  to  keep  the  temperature  below  600K  in  order  not  to 

destroy  the  sample  which  visibly  darkened  by  only  briefly  heating  it 

lie 

above  that  limit.  Therefore  we  can  assume  <<  1  *  and  Ed-  1~B 

reduces  to 

-  8irhcX-,exp(-  (2-B) 

The  total  intensity  I  in  the  spectral  range  of  sensitivity  of  the 
photomultiplier  tube  is  given  by 
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A  2 

Ib  a  /  K(A)UAdX  (3-B) 

1  i 

where  K ( A )  Is  the  quantum  efficiency.  It  turned  out  to  be  a  linearly 
decreasing  function  of  the  wavelength  between  600nm  and  650nm,  the 
wavelength  interval  utilized  in  our  experiments.  This  yields 

I  a  exp(-  (4-B) 

AkT 

where  A  is  at  the  center  of  this  sensitivity  region  which  is  determined 
by  the  photomultiplier  and  a  short  wavelength  cut-off  filter  (see 
below ) . 

In  order  to  eliminate  the  unknown  proportionality  constant  in  Eq. 
4-3,  one  can  utilize  the  ratio  of  the  first  derivative  with  respect  to 
time  and  this  expression.  This  ratio  i3  readily  determined 
experimentally: 


(|i)/I  a  (a/T* )exp(-a/T) (dT/dt ) 


(5-B) 


and 


(^1)/I  -  (a/T£)( dT/dt) 


(6-3) 


where  a  -  — . 

Ak 

The  time  dependence  of  the  temperature  is  given  by  Eq.  8  of 
Section  III. 2.  We  may  consider  the  simplest  case  which  for  the 
temperature  is  proportional  to  the  time  t.  This  is  valid  only  for 
4at  <  w2  and  yields 


T-T  a  Bt/w2 
o 


(7-B) 


I 


I 
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Then 


(ii)/I  -  [a^Bt/w2  ♦  Tq)2](B/w2).  (8-B) 

The  quantity  B  is  calculated  by  determining  this  ratio  experimentally 
and  careful  measurement  of  the  1/e  full  width,  2w,  of  the  laser  beam. 

EXPERIMENT  RESULTS 

One  of  the  4  different  dosemeter  conf igurations  wa3  a  polyimide 

film  loaded  with  LiF  powder.  Both  materials  of  thi3  84ym  thick 

composite  film  absorb  the  10. 6pm  heating  photons  only  moderately.  One 

can,  therefore,  assume  that  Eq.  8  is  applicable. 

The  parameter  combination  B  of  Eq.  9  is  straightforwardly 

determined  by  recording  the  time  evolution  of  the  incandescent  emission 

at  a  narrow  wavelength  interval  around  A.  Using  a  step  filter  having  a 

low-wavelength  transmission  cut-off  at  600nm  (interference  filter  type 

6135-650,  Optics  Technology,  Palo  Alto,  California).  Together  with  an 

EMI  9924B  PMT  yielded  A  -  625nra.  The  incandescent  signal  obtained  by 

heating  this  foil  with  P  -  4.8  Watts  and  2  *  0.11 8cm  is  3hown  in 

Fig.  1-B  together  with  its  derivative  as  obtained  by  digital  processing 

with  the  Data  Pecision  DP60QQ  waveform  analyzer.  From  the  ratio 

(■*i)/I,  again  calculated  directly  with  the  DP6000,  B  wa3  determined  to 

be  27cm2  K/s  for  the  laser  power  stated  above.  The  condition  4at<w2  is 

fulfilled  with  the  values  a  -  0.0175cm2/3  (see  Section  III. 2),  t  »  .03 
0.0  7  ^ 

sec  and  w  -  Q^jbcHC. 

It  should  be  pointed  out  that  the  above  method  is  not  applicable 
for  the  dosemeter  configuration  described  in  Sections  III .  1 ,  3,  and  4 
because  here  the  z-dependence  of  the  temperature  rise  cannot  be 
neglected  and  Eq.  8  i3  not  valid. 
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APPENDIX  C 


AN  ADVANCED  RESEARCH  FACILITY  FOR  LASER  HEATING  IN  THERMOLUMINESCENCE 

After  completion  of  the  research  on  the  heat  transfer  from  a  Gaussian 
laser  beam  to  semi-infinite  thin-layer  dosimeter  configurations  and  of 
the  resulting  thermolurainescence  response  curves,  a  new  apparatus  was 
designed.  This  facility  is  needed  to  study  these  aspects  in  a  situation 
that  closely  resembles  the  experimental  conditions  which  have  evolved  as 
the  ones  that  most  likely  will  exist  in  actual  field  use  of  laser-heated 
thermoluminescnece  dosimetry  readers  to  be  developed  in  the  next  two 
years. 

These  conditions  are  the  following: 

a.  The  laser  beam  must  have  a  uniform  intensity  profile,  be  highly 
stable,  reproducible,  and  computer-controllable. 

b.  The  laser  power  must  be  at  least  40  Watts. 

c.  For  maximum  sensitivity  of  the  dose  measurements,  the  spot  size  should 
be  at  least  3x3  mmi  The  high  laser  power  is  needed  to  heat  this  spot  in 
less  than  200  ms  to  400  deg.  C. 

d.  Since  the  beam  optics  exhibits  a  polarization  dependence  of  the 
reflectivity,  it  is  mandatory  to  assure  that  the  beam  polarization 
remains  fixed  during  operation  of  the  system. 


The  new  system  that  was  constructed  and  is  presently  being  tested 
features  a  40  Watt  RF  excited  CO  laser  emitting  a  one  cm  diameter 
"halo"  beam  (doughnut-shaped  cross  section)  of  fixed  polarization  (Model 
L35PS,  Directed  Energy,  Inc.,  Irvine,  California).  The  polarization 
stability  is  achieved  by  an  intra-cavity  ZnSe  polarizer  plate  that  is 
mounted  inside  the  beryllia  discharge  tube.  The  uniform  beam  profile  is 
obtained  by  tightly  focussing  this  beam  with  a  short  focal  length  Ge  lens 
into  a  gold-coated  optical  channel  of  3x3  mm’cross  section.  The  optical 
layout  is  schematically  shown  in  Fig. 17. 

A  ZnSe  beamsplitter  is  placed  immediately  after  the  focussing  cube  (shown 
as  a  concave  mirror  in  Fig. 17),  and  its  position  is  adjusted  so  as  to 
direct  the  small  reflected  beam  to  the  center  of  the  PZT  detector.  The 
chopping  wheel,  turning  at  2000  rpm,  has  three  1  mm  wide  slots  that 
permit  the  beam  to  hit  the  detector  once  every  10  ms.  The  PZT  detector 
has  a  5x5  mm’sensitive  area.  The  pulsed  signal  it  generates  is  used  as 
the  input  signal  for  the  laser  control  electronics  described  below. 

A  nitrogen  purge  inlet  is  located  just  below  the  small  Ge  lens  at  the  top 
of  the  optical  channel,  and  the  nitrogen  flows  down  onto  the  dosemeter 
and  into  the  dosemeter  reading  chamber.  The  dosemeter  is  positioned  at 
the  exit  of  and  very  close  to  the  optical  channel.  A  light  guide 
transmits  the  emitted  luminescence  to  the  photomultiplier  tube  (PMT). 


The  Control  Electronics: 

When  a  dosimeter  read  cycle  is  initiated,  the  following  operation 
sequence  is  performed  by  the  control  electronics 

1.  The  laser  is  turned  on  and  the  laser  power  is  measured  and  adjusted 
until  the  laser  is  stabilized  at  the  desired  polarization  and  power 
level.  This  typically  requires  100  to  300  ms.  During  this  time  the 
shutter  remains  closed  and  the  dosemeter  is  not  exposed  to  the  laser  beam 

2.  When  the  desired  laser  power  is  obtained,  the  shutter  is  opened  for 
the  desired  heating  time.  It  is  possible  that  the  laser  power  and/or  the 
polarization  may  suddenly  change  during  the  time  that  the  dosemeter  is 
heated.  If  this  occurs,  indicators  on  the  front  panel  of  the  control 
electronics  are  set  to  indicate  the  occurrence  of  this  condition.  A 
trigger  output  provides  an  output  pulse  corresponding  to  the  heating 
cycle 

3.  At  the  conclusion  of  the  heating  cycle  the  shutter  is  closed  and  the 
laser  is  turned  off.  The  system  is  immediately  ready  to  begin  another 
read  cycle. 

Laser  power  control  is  accomplished  by  pulse  width  modulation  of  the 
30  kHz  laser  modulation  signal.  The  laser  power  control  function  will 
normally  maintain  the  measured  power  within  less  than  1%  of  the  laser 
power  setpoint.  However,  a  sudden  change  in  the  laser  polarization  or 
output  power  may  result  in  a  deviation  of  more  than  1%.  To  check  for  the 
occurrence  of  this  condition,  a  power  checking  circuit  is  used  to  check 
for  a  power  error  of  greater  than  3%.  The  dosemeter  read  cycle  is 
inhibited  when  this  condition  occurs.  If  a  dosemeter  read  cycle  has 
already  been  started,  an  indicator  is  set  to  warn  of  the  power  error. 

In  addition  to  the  power  checking,  there  is  another  checking  circuit 
that  is  used  to  avoid  the  use  of  the  laser  when  undesirable  polarization 
is  present.  Should  the  beam  polarization  deviate  from  the  desired 
polarization,  the  beamsplitter  reflects  a  greater  or  smaller  percentage 
of  the  incident  power.  This  gives  an  erroneous  indication  of  increased 
laser  power.  This  condition  can  be  detected  because  it  results  in  a 
substantial  change  in  the  control  voltage.  For  a  given  laser  power 
setpoint,  an  acceptable  range  in  control  voltage  can  be  determined,  and 
these  limits  are  then  used  to  check  for  this  polarization  change.  This 
checking  function  is  referred  to  as  control  voltage  checking.  A  dosemeter 
read  cycle  is  inhibited  when  this  condition  is  present,  and  a  control 
error  indicator  is  set  if  this  condition  occurs  after  a  dosemeter  read 
cycle  has  been  started. 

The  addition  of  a  temperature  stabilized  base  for  the  present  laser 
head  (Neslab  refrigerated  recirculating  heat  exchanger  model  RTE-5DD) 
further  assures  stable  operation  of  the  laser  as  it  maintains  the  laser 
discharge  tube  at  a  constant  temperature  (±.1  deg.  C). 


Optica]  chamber 


fig.  17:  Optics  of  the  Proposed  Laser  TS.I)  Reader  (Vertical  layout) 
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